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Abstract 

We study the behavior of solutions to a Schrodinger equation with large, rapidly oscil- 
lating, mean zero, random potential with Gaussian distribution. We show that in high 
dimension d > m, where m is the order of the spatial pseudo-differential operator in 
the Schrodinger equation (with m = 2 for the standard Laplace operator), the solution 
converges in the L 2 sense uniformly in time over finite intervals to the solution of a de- 
terministic Schrodinger equation as the correlation length e tends to 0. This generalizes 
to long times the convergence results obtained for short times and for the heat equation 
in [2]. The result is based on the decomposition of multiple scattering contributions 
introduced in [6]. In dimension d < m, the random solution converges to the solution 
of a stochastic partial differential equation; see [1, 13]. 

1 Introduction 

There is a long list of derivations of macroscopic models for solutions to equations involving 
small scale heterogeneities £<1. One very successful framework is that of homogenization 
theory. In the limit e — > 0, it is shown that the heterogeneous solution converges to the 
deterministic solution of an effective medium equation [5, 9]. There are cases, however, 
where the solutions in the limit e — > remain stochastic, typically in low spatial dimensions; 
see e.g. [1, 7, 8, 11]. 

The results of [1, 11] apply to parabolic equations of the form of a heat equation with 
large, mean-zero, highly oscillatory, random potentials (zeroth-order terms). In [2], it was 
shown that for large dimensions, the random solution converged to a deterministic solution, 
which is consistent with the homogenization framework. However, such results could only 
be obtained for short times, and it is unclear whether they hold for larger times. 

In this paper, we revisit the homogenization limit for the (time-dependent) Schrodinger 
equation. Because the solution operator is unitary in this case, we expect to be able to 
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control the long-time asymptotic behavior of the solution. The method of proof, as in [1, 2], 
is based on a Duhamel expansion of the random solution in terms involving increasing 
numbers of scattering events. As the number of terms grows exponentially with the number 
of scattering events, we need to assume that the potential is Gaussian in order to control 
such a growth. The summation technique used in [2] cannot extend to long time controls. 
A similar difficulty occurs in the derivation of radiative transfer equations for the energy 
density of high frequency waves propagating in highly oscillatory media [4, 6, 12]. A precise 
summation of the scattering terms was introduced in [6] to allow for long-time expansions. 
We adapt this technique to the asymptotic analysis of Schrodinger equations with large 
potentials. 

We now present in more detail the model considered in this paper and the main con- 
vergence result. Let m > 2. We consider the following Schrodinger equation in dimension 
d > m: 

igi + (P(D) ~ -^q^)))u e (t,x) = 0, t>0,xGR d ^ 

u £ (0,x) = u Q (x), x G R d . 

Here, P(D) is the pseudo-differential operator with symbol j5(£) = |£| m . We assume that 
q(x) is a real valued mean zero stationary Gaussian process defined on a probability space 
(fijj 7 , P), with correlation function R(x) = K{q(y)q(x + y)}, and the non-negative power 
spectrum is radially symmetric, smooth, and decays fast. For simplicity, we assume 
R G S(]R d ). In fact, ||-R||i2d,i2d < +oo, where 

ll/lk* := IK*>*<V.>*/(*)|| a , <*) := (1 + x 2 ) 1 ^, (1.2) 

is enough. We choose the initial condition uq(x) to be smooth such that uo(£)(£) M ^ L 2 (R d ). 
For any finite time T > 0, the existence of a weak solution u £ (t, x) G L 2 (S1 x R d ) uniformly 
in time t G (0, T) and < e < £o can be proved by using a method based on Duhamel 
expansion. 

As e — > 0, we show that the solution u £ (t) to (1.1) converges strongly in L 2 (QxR d ) uniformly 
in t G (0, T) to its limit u(t) solution of the following homogenized equation 

(i^ + P(D)-p)u(t,x) = 0, t>0, xeR d 



where the potential is given by 



u(0, x) = uq{x), x G 



/ 



dC (1.4) 



lei 

The main result of this paper is the following convergence result: 

Theorem 1.1. There exists a solution to (1.1) u e (t) uniformly in < e < e$ for t > 0. 

Moreover, we have the convergence results for all t G (0, T) 

limE||(« e -«)(*)||l = 0. (1.5) 
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The rest of the paper is organized as follows. Section 2 recasts the solution to (1.1) as a 
Duhamel series expansion in the frequency domain. We estimate the 1? norm of the first no 
terms by calculating the contributions of graphs in three categories similar to those defined 
in [2]. Section 3 estimates the L 2 norm of the error term by first subdividing the time 
integration into time intervals of smaller sizes, and then using Duhamel formula in each 
time interval. This method, introduced in [6], significantly improves the error estimates 
compared to the direct estimates of infinite Duhamel terms and enables the elimination of 
the restriction to short times. The estimates given in these sections are used in section 
4 to characterize the limit of the solution u e (t,x). Section 5 provides the proofs for the 
inequalities used for justifying the estimates in the previous sections. 

The analysis of a parabolic equation of the type of the heat equation (with idt replaced 
by dt) is performed in [2] for d > m. Up to a logarithmic correction, we expect the limit of 
the solution to (1.1) to be deterministic also for the critical dimension m = d. In [2], the 
random fluctuations about the deterministic limit are also analyzed for the heat equation. 
We expect a similar behavior to occur for the Schrodinger equation (1.1) for short times. 
We do not know the behavior of the random fluctuations for arbitrary times t £ (0, T). 

In lower spatial dimension d < m, the limit of the solutions to (1.1) as e — > remains 
stochastic. This behavior was analyzed for the heat equation in [1, 10]. The limit of u £ is 
then shown to be the solution of a stochastic partial differential equation with multiplicative 
noise (written as a Stratonovich product). The analysis of (1.1) for d < m is performed 
in [13]. Note that several results of convergence may be extended to the case of random 
potential with long range correlations or random potentials that display both temporal and 
spatial fluctuations [3]. 

2 Duhamel expansion 

We denote by e ltH the propagator for the equation (1.1). The Duhamel expansion then 
states that for any no > 1 



no — i 




(2.1) 



n=0 



where for H := (— A) 2 , we have defined 



u n ,s) : = Hn-^rl 



■t 



, t n n 



(J] ds kW - £ s k )e is ° H °q(-) ■ ■ ■ q(-)e 



iSnHo 



n , (2.2) 



fc=0 fc=0 



and the error term is given by 




(2.3) 



We shall choose 



n = n (e) : 



7l loge 



(2.4) 



log I log £ 



3 



for some fixed < 7 <C A sufficiently small, where A is defined as 

. f d - m m < d < 2m , n 
{ m d > 2m. (2 - 5) 

For any subset / C N, we define the kernel for the evolution in the Fourier space as 

Jo kei kei fee/ 

Hereafter, we use the notation £ m = |£| m . In the special case where I = {0, • • • , n}, we 
denote by K(t;£,n) := K(t;£,{0,- ■ ■ ,n}) and £ n := £ Jn . Denote £ n6 = £/ n \{ }- 
Let us introduce q e (Q = e d ~^q(e^), the Fourier transform of e - '2q , (|). Denote the contri- 
bution from the potential term by 

n 

L&n) :=!!&(&-&-!). (2.7) 

fc=i 

We may rewrite the n th order wave function as 

« e ,n(t,&) = J K(t;Z,n)L(£,n)u (Z n )d£ n (j . (2.8) 
We need to introduce the following moments 

t/ £ "(«o) = ER»}> (2-9) 

which are given by 

t£(Uo) = {~i) n f K(t,Z,n)E{L(S,n)}u (tin)dZ n>i i, (2.10) 

and 

[/ e fl - m (Uo,Co)=lER, n (t 1 0^(i,C)}, (2.11) 

which are given by 

U?' m (t,£o,(o) =(-i) n+m J K(t, C, n)K(t, C, m)E{L(t n)Z(^}fi (^)^(U)d€ n) orfC m ,6- 

(2.12) 

We need to estimate moments of the Gaussian process q e . The expectation in Jj^' m 
vanishes unless there is n S N such that n + m = 2n is even. The moments are thus given 
as a sum of products of the expectation of pairs of terms q e (£,k — £fc+i)i where the sum runs 
over all possible pairings. We define the pair (£&,£/), 1 < k < I, as the contribution in the 
product given by 

- 6)4(6-1 - 6)} = e d - m %(£fe - - ffc-i + 6 - 6-1) 

= e d - m r(e(4 - &_i))r(e(& - - + 6 - 6-1) 

(2.13) 
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with r(£) := R(£) 1/2 - 
Define 

n 

F(€, n) := ]J r(e(£ fc - &_i))«o(£ n )- (2-14) 
fc=l 

Denote by A„- the product of delta functions associated with the pairing ir. Our analysis is 
based on the estimate of 

7r6ll(n) 

with 



I* ■= {~i) n J K(t,tn)A w (£)F(£,n)dt (2.16) 

Ju?> m (t,a,odn= Y, c - ( 2 - 17 ) 



and 

7r£il(n,m) 

with 

Cn : = (-*)™+ m y K(t, t n)K(t,C,m)A 7T (t C)£(£o - Co)F(£, n)F(C,m)^ ft dC n . (2.18) 
By Lemma 5.1, i^(t; £, /) can also be written as 

K(t; /) = fe** / doe"** J] * . (2.19) 

We let r\ = t~ 1 in this section. Therefore, I n in (2.16) and C n in (2.18) can be written 
explicitly as 

I* = / doe"*- f[ ! ^fe%^^^)^(e.)^ > 6, (2-20) 

J , . "TUT"/ 



and 



n ^ m .. n 

Cw = -in+m e 2t V / dad p e -i(«-P)t JJ JJ J] r ( e (£fc - fc-l)) 

^ fe = o « + & + "7 t 13 + C ™ + 111 k=i 

m 

Jr(e(Ci - 0-i)A^,C)u (£n)£(C^(£o - Co)d&»dCm- 



x 

=i 



(2-21) 

In order to consider the two sets of momenta in a unified way, we introduce the notation 

_ J a < k < n , , 

afc ~\/3 n + l<fc<n + ro + l, 1 j 

and define = Cm-k f° r < /c < m. (2.21) can then be rewritten as 

/n+m+l -. n+m+1 

jti ^ + ^ + ^^4^ (2.23) 

^o(£n)^o(£m)<5(£o - £n+m+l)d£ re+m+1 . 



In each instance of the pairings, we have ft terms k and n terms I = l(k). Note that 
l(k) > k+1. The collection of pairs (£k->£l(k)) f° r ^ values of k and n values of l(k) constitutes 
a graph ir E II(n, m), where n(n, m) denotes the set of all graphs tt with n copies of q and 
m copies of q. The graphs are defined similarly in the calculation of E/"(i,£o) iu (2.10) 
for n = 2n, and we denote by Il(n) the set of graphs with n copies of q. We denote by 
Aq = Aq(tt) the collection of the n values of k and by Bq = Bq(tt) the collection of the n 
values of l(k). 

Now we introduce several classes of graphs for C n . We say that the graph has a crossing 
if there is a k < n such that l(k) > n + 2. We denote by H c (n,m) C U(n,m) the set of 
graphs with at least one crossing and by H nc (n,m) = H(n,m)\Il c (n, m) the non-crossing 
graphs. We denote by simple pairs the pairs such that l(k) = k + 1, which thus involve a 
delta function of the form — Cfc-i)- The unique graph with only simple pairs is called 

the simple graph, which is denoted by H s (n,m). H ncs (n,m) = H nc (n,m)\H s (n,m) denotes 
the set of non-crossing, non-simple graphs. 

We also use the notation IL s (n) and n ncs (n) defined for which denote the simple 
graph and the set of non-crossing, non-simple graphs, respectively. 

We shall estimate F(£,n) before we proceed to analyze the graphs. For the initial 
condition, we define 

HO = <0 M M0I- (2-24) 
From our assumption on no given in Section 1, we have $(£) £ L 2 (R. d ). Since ||-R||i2<z I2d < 
+oo, we have 

C_ 

Hence, we obtain the estimate for F(£,n) 



K£) < tJvm- ( 2 - 25 ) 



|^,n)|<n 



1 m« 



k= 
n 



< n 1 to 



k=l 



<e»> M ti(^> M ii<e(&-&-i)> M ' 

We have the freedom to choose /i and Z2 between and n. 

Analysis of the crossing graphs 
Lemma 2.1. If tt £ TL c (n,m), we have 

|CV| < (Cloge)V||$(e)||l. (2.27) 

Proof. Denote by (£g m ,6(g m ))) 1 < ?n < M, the crossing pairs and define Q = m&x m {q m }. 
Let us define A' = A Q \{Q}. From (2.26), we have 
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The terms 



|Q, fc+ gm +it -i| for k £ A' U {0, n, n + 1, n + m + 1} are bounded by C. This allows 



us to obtain 

n 



i 



i 



i 



a + ^m + i t -l| |0 + £« + j t -l| ( eeo )2d 

— 2^(6: ~ + &(fc) ~ £z(fc)-l) 



<5(£q - £q-1 + £j(Q) - &(Q)-l) 
| a + e + , t -l M/3 + e + ^l|l^)l^ 



(2.29) 



n+m' 



For each iGi'U {0}, we perform the change of variables — )• — , and define 



% kA'u{0}. 



(2.30) 



We then find the estimate 

\C n \ <C n J dad(3 

n 



,2d 



<*( — - £1-1 + - £f(fc)-lM£Q - CQ-l + - &(Q)-l) 

^T2ll^(^)| 2 ^ n+m . 



(2.31) 



We now estimate the above product. Assume Q < n and n = l(ko). When Q = n or 
KQ) = n + m + 1, the derivation of the same estimates is simpler and left to the reader. 
Define k\ such that l(ki) = n + m + 1. For each k S A'(7r)\(ko U fci), we use (5.25) below 
to find the estimate 

1 1 £ 

(2.32) 

The integration in the £z(Q) variable is estimated by using the above delta function. The 

delta function for k = ko G ^'(tt) may be written in the form <5(£q — + £o — + 
iAf-1 



]Cm=i £? m ~~ £?m-l)> an d is thus used to integrate in the variable £q. The term 
is used to integrate in the variable £fc . The integral in a and j3 is estimated using (5.3) by 



-da < C-, 



The integral in £o is estimated using (5.18) by 



1 + log. 



1 


1 + log+ 




- 2 


<£o} M 


ie m /e m - 


e + 


2 



< Ce A (e^ v l). 



(2.33) 



(2.34) 
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Following the usual convention, we use a V b := max{a, 6} and a A b := min{a, b}. The delta 
function <5(£ fcl - - Cn+m+i - £n+m) is seen to be equivalent to 5(£ n - £ n +i), which 

handles the integration in the variable £ n +i- Finally wg integrate in ^ n 



I 



V l)$(en)| 2 ^n < II^COIll, (2-35) 

and obtain 

|C w |<(C|loge|)V||*(0||!. (2-36) 

□ 

Using Stirling's formula, we find that |IT c (n,m)| < 2 ^{n-iy. ^ s bounded by (^r) n . After 
summation in n, m < uq, we obtain 

n n — lno— 1 r>- 

E||$> n , £ -^)(i)||I<£ ^(^f(C|loge|f £ A ||$(e)||| 

n=0 n=0 m=0 

<n^(C|log e |re A ||$(e)||i< £ A - 3 ^, 

where a < 6 means a < Cb for some C > 0. 
Analysis of the non-crossing graphs 

Lemma 2.2. If it £ H ncs (n), we have 

{I^KiClogs^e^mOl (2-38) 
where 0<iJ€l. Moreover, if n & H ncs (n,m), 

\C W \ < (Clogefe 2X ^\MO\\l (2-39) 
Proof. In a graph ir 6 II ncs (n), the delta function 

5(£o " Cn) (2.40) 

is obtained by adding up all the delta functions in A^. We perform the change of variables 
for all k £ Aq, ^ — > — and define as before 



,e _ I k i Aq 

% keA . (2 - 41) 



We shall solve the following two cases in different ways. 

(i) If there exists ki € Aq, such that for all k satisfying ki + 1 < k < l{k2) — 2, k £ Aq, 
are simple pairs, then the delta function 8(-^ — Cz(fe 2 )-i) * s present. From (2.26), 

we have 

F( M < C"WWI^ II (^p- C™> 
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The estimate of integration in £fc 2 is then obtained by using (5.18) below: 

The delta function in which is involved is equivalent to <5(£ n — £o) ? which we use to 
integrate in £ n : 

For k S Aq , k ^ k2 , we have the estimate 

/ -,g l) *| tt + g /]. + <t V <7 " + f '"> - Ww&m s ci iog4 

(2.45) 

The estimate of integration in the variable a is then given by 



Mi-*) 
am vi 



— <C$ L - S) . (2.46) 



The extra term ^ that arises in the last estimate can be canceled by the term l/(£o) 2<i 
in (2.44), which concludes (2.38). 

(ii) If there exists no such &2 £ satisfying the condition in case (i), then we first 
delete all simple pairs that exist in the graph ir. In fact, the simple pairs can be handled 
first by using the bound as in (2.58). Therefore without loss of generality, we need only to 
consider a graph n with no simple pair. Let us define k± = min{fc|A; € Ao,l(k) — 1 £ ^4o}> 
and k$ = l(k$) — 1. Note that k$ > fej + 1. 

We have from (2.26) that 
The integration in £j(fc 4 ) provides the terms which we will need for integration in ^ 5 



We can estimate the integration in ^ 5 using (5.18) 

„ 1 1 1 



(2.49) 



|a + |6 5 -6 4 +ea 4 -il m /e m + ^ 
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All the other integrations are handled the same way as in case (i). In order to make sure 
that no integration above is affected by other integrands we plan to use for integrating in 
other variables, we just need to first integrate in with index in decreasing order and 
then integrate in with index in decreasing order. This gives (2.38). 
If 7r G II ncs (n, m), we may denote the pairings for k < n and for n + l<A;<n + m + lby 
7Ti and 7T2 and then ir = tti U tt2, since there is no crossing in ir. Hence it follows that 

lei < J ii^o^m < (ciog^v^p^m. (2.50) 

□ 



Similarly to (2.37), we obtain 

no — I no— lno — 1 „_ 

ii ( u ? - Ksmwl < E E(-) n ( <7 i 1 °g £ D fie2A(1 ^ ) ii $ ^)iii 

n=0 n=0 m=0 6 (2.51) 

<n^(C|loge|) no e 2A(1 - 5) ||$(OHi < e 2A(1 ~ 5) ~ 37 , 

where 

u% ■= i« ( 2 - 52 ) 

for 7r G H s (n). 

Collecting the results obtained in (2.37) and (2.51), we have shown that 

n -l 

E|| E (ue,n - UZ.)(t)\\2 Z e 2X ^~ 3 \ (2.53) 

n=0 

Analysis of the simple Graphs 
Lemma 2.3. If tr G Il s (n), we have 

/-in 

%\ < (j^ + O(C n e m / 2 ))\u (a (2-54) 
(n/2)\ 

Moreover, if tv G U s (n,m), we have 

/-m+m 

™ £ V/2)K,n/2)! + °( C " +m£m/2 ))ll»»«)H- ^ 

Proof. In the case of a simple graph tt, we can explicitly write out the product of delta 
functions 

A. = J] - (2.56) 

k€A 

which is independent of for all k G Aq, and forces = £o for an & ^ Re- 
integrating in £fc for all k G Bo using delta functions, we obtain 

I n = f K(t,Co,--- ,to,Hk ai ,--- ,Ck an/2 ) I] ^' m R(^k-^o))M^ Ao , (2.57) 
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where {k ai , • • • , K n/2 } = A . 
This implies 



g ^R( £ (g fc -g )) 



Define 



K(£o)|- (2.58) 



«,^o) = / ■ ml- ( 2 - 59 ) 

Perform the change of variable — >■ ™. This gives 

0ar,(£o) = I R ^ k ~ ^ dfo. (2.60) 
It is clear from Lemma 5.8 that 

|e a ,„(£o)| < C. (2.61) 

Thus (2.58) already implies 

\U\ < C n \M0\- (2-62) 

However, this estimate is not sufficient and there is in fact a term l/(n/2)! missing. We 
now recover this factor. 



Introduce the notation 



9(0 = Hm 0^(0- (2.63) 



From the estimate in (5.36), we have 

|0a lJ? (eo) - e(&)| < e m/2 (|a + e m |kr 1/2 + H 1/2 )- (2.64) 
We shall show that the leading term of I n is 

K(t,£o,--- ,eo)e(^o)" /2 no(eo). (2.65) 
In fact, the error term is bounded by 



da- 



h m/2+1 - e n/2 m\MZo)\- (2.66) 



a + £o + iv\ n/2+ 



Using the uniform bound on in (2.61) and the estimate in (2.64), we can bound (2.66) 
by 0(£ m / 2 C"|n (£o)|). 

We can now use Lemma 5.1 to bound (2.65) by 

«o(£o)|- (2-67) 



(n/2) 



Finally, we obtain (2.55) as an immediate consequence of (2.54). This concludes Lemma 
2.3. □ 
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3 Partial Time Integration 

In this section, we estimate the I? norm of ^ no ,e- The central idea is to subdivide the time 
integration into smaller time intervals of size t/nie) with 

K{e):=\\oge\ lh2 . (3.1) 

We then use the Duhamel formula to estimate the evolution in each time interval. 

-no 



Recall the error term ^n ,e = ^2n=n u n,e- The Duhamel formula states that 



*GW(t) = (-1)4/2 f d.se^ H q(-)u no ^ £ {s), (3-2) 



lo 

where e ltH denotes the propagator of equation (1.1) 
Let Oj = jt/n for j = 0, 1, ■ ■ ■ , k. Rewrite 



= H)~^ Y,^' 6 ^ I ^ e ^ H q{-)u no . h£ { S )ds. (3.3) 

3=0 



Define the n-th term of the Duhamel expansion for the operator e^ j+1 S ^ H in (3.3) as 

1 r e j+i ft ft n ~ no n ~ n ° 

J9< Jo Jo tt) (3.4) 



e isoH °q(-) ■ ■ ■ q(-)e is —o + ^o Uno _ l £{s)ds _ 

We may further obtain the form of w n ,no,0j in terms of uq by writing ii no _i i£ (s) out explicitly 
using (2.2) 



i r e j+1 r t r t n 
s^' J 6j Jo Jo 

n— no n 



=H) n (^r/ " / ■ / (11**) 

6j JO JO k=Q 



(3.5) 



6(6 j+1 -s-J2 s k )S(s - ]T Sk y s ° H °q{- £ ) ■ ■ ■ q(-)e^ H °u ds. 

k=0 k=n— no+1 



Uin^nofijix) = -^q(-)u 4 no-l,n ,e 3 (x). (3.6) 



The amputated versions of these functions are defined as 

1 ,x s 

Then Duhamel formula then gives 

1W = Ui + U 2 , (3.7) 

where 



K-l 

no<rt<4no 3=0 

U 2 (t) = (-i)J2e^- e ^ H f 3+1 u ino , no , ej (s)ds. 



(3.8) 



3=0 
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From the unitarity of e^* d i+ 1 ) H and the triangle inequality, we can bound U\ by 

ll^llli < ChqK ^2 \\ u n,no,Sj(^3+l)\\l ^ Cn 2 n 2 sup H'Wr^nojfj \\l (3-9) 
- . n nn<n<4nn 

Applying the Cauchy-Schwarz inequality, we can bound U2 by 

\\U2\\\<t2_\ I \\Uin ,n ,eAs)\\lds < t 2 SUp SUp ||^4n ,no,% ( S ) 111- ( 3 - 10 ) 

Denote 

J n _ no+ i jn := {n - n + 1, • • • , n}. (3.11) 

Define the free evolution operator with constraint given by the parameters no and 9 as 

K*(9j + i,0j;£,n, n ) := / K(9 j+1 - s; £, n - n )^(s, £, / n -n +i,n)ds- (3.12) 

We can write the wave function in Fourier space Un,n ,0j as 

,n,n )L(€,n)u (€n)d€ n Q. (3.13) 

We then write 

E||u w ,(0 i+1 )|| 2 = J] C# (3.14) 

7rSlI(n,n) 

where 

C#:= y d£dCK#(p j+1 , 0f,£, n, n )K#(9 j+1 , 9 f , (, n, t*>) A ff F(£, n)F%n). (3.15) 
For the amputated function, we have 

E||u w ,(0 i+ i)|| 2 = J2 d *> ( 3 - 16 ) 

7rGlI(n,n) 

C*:= J d£dtk#(0 j+1 , Ofit, n, n )K#(9 j+1 ,9 j ; C, n, n )A 7T F^, n)Ffcn), (3.17) 
where 

K*(6 j+ i,6j;£,n, n ) := / - s; £, Ii^-^Jif (s, £, 7 n _ no+ i in )ds. (3.18) 
Recall Lemma 5.3. We can extend it to the following identity for K#(9j + i,9j-,£,n,no): 

r8j +1 r+00 

K*(9 j+1 ,9 j -^,n,n ) =- / dse^'^e^ \ dadae'^ 8 ^ 1 - 3 ' ] e~ i&s 

J6j J-00 

n-no ^ ^ (3.19) 

X 



P a +^+^ fe6 /^ +1 «+^ m +^r 

<J ""t 1 71 — 7lQ-|-l,7l 
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where we choose r\ := (i/ft) \fj := t . We can integrate in s to have 

e -i8 j+1 a+9 j+1 ij _ e -ia(9j +1 -6j)-i9ja+(6j +1 -6j)ri+6jfj 



/+oo 
dada- 
-oo 
n— no 

X 



n 



a + £{" + ir, 
k=o ^ k ' kel 

Hence we can bound C^f by 

/r+oo 
di / dadadfidp 
J — oo 

1 



n 



a — a + i(r] — fj) 
1 



n — tiq + l,ra 



n 



0<fc<2n+l 1 ' ^ k ' k=l,k^n+l 



a - a + i(rj - fj)\ \p - p + ifr - 

2n+l 

r(e(&-&_i)), 



where and r]k are defined as 



OA: 



Vk 



'a if k < n — uq 

a if n — uq + 1 < k < n 

P ifn+l<k<n + n,Q 

P if k>n + n + 2 

fj if k < n — no 

i] if n — n,Q + 1 < k < n 

i] if n + 1 < k < n + no 

fj if k > n + no + 2. 



We present the following lemmas to show limE||\I/ no e || 2 

e— >o ' 



Lemma 3.1. Let n = 4tlq. For any ir £ U(n,n) we have 

(C|loge|) 4no 



0. 



\c*\ < 



K n 



Proof. The following bound can be easily obtained by using Lemma 

|C#l<(C|loge|) 4 'T 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



To recover the denominator in (3.24), notice that among the % for k S Bq, there are at 
least n — 2riQ — 2(> no) of them with rjk = n/t. Hence 



n 



i 



< 



< 



K n 



(3.26) 
□ 
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Lemma 3.2. If ir E II c (n, n), then we have 



\C*\ < {C\ogeTe x \\m\\ 



2- 



(3.27) 



Proof. The proof of Lemma 3.2 is essentially the same as in Lemma 2.1. The only difference 
is that the integral in a and a (/3 and f3) is estimated using Proposition 5.3 by 



1 



C 



\a — a + i(r/ — fj)\ \a + £™ + w?| |a + ^/e 1 " + ir? 
and the integral in £o is estimated using Proposition 5.6 

1 + log+ 



dada < 



1 + log H 



is? - e m /e m + "j| 



1 



tin /_m tin 



-d£o < Cs x (& V 1 



<6> 2d \^/e m -^ + iv\ 2 



Lemma 3.3. If tt E n ncs (n, n), £/ien we /i<roe 



|cf | < (ciog^^^-^H^coil 



2- 



(3.28) 

(3.29) 
□ 

(3.30) 



Proof. The proof of Lemma 3.3 is similar to that of Lemma 2.2. The only difference is that 
the integral in a and a (/3 and f3) is estimated by 



A(l-5) 

(a~ — V 1) 



<C / (a^~ V 1) 



1 



1 



1 



\a — a + i(rj — fj)\ \a + £Jf + Z7y| |q + ^ + 277] 



dada 



\a + $ + ifj\ 



log H 



+ 1 



A(l-5) 



Lemma 3.4. 7/7T E n s (n, n), i/ien we Ziaue 



(3.31) 



□ 



(3.32) 



Proof. The proof is essentially the same as in Lemma 2.3. We shall not repeat the argument 
here. □ 

We now apply the above lemmas to estimate *$> no e . From Lemma 3.2, 3.3, and 3.4 we have 



From Lemma 3.1 we have 



K n ° e 



(3.33) 



(3.34) 
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The L 2 estimate of ^ no ,e is therefore given by 

oo 

E|| J2^At)\\l = n^noJl < 2(E||C7 1 ||1 + E||C/ 2 ||1) -+ (3.35) 

as e — > 0. 

4 Homogenization 

We come back to the analysis of U £)S {t, £). We find that U EjS is the solution to the following 
equation 

U £ , s = e u ^n (0- fe' sr f ' e iS ^ / e d - m R(e(Ci - t))U e ,.(t - s - ai , O^idsds, 
Jo Jo J 

(4.1) 



Lemma 4.1. Let us define il e (i,£) to be the solution to 



with p e = j Rd R ^ m dt;i . We have the convergence 

W^-iXeXt)] <max{e m ,e d - m |loge|}. (4.3) 
Proof. (1) We obtain from Duhamel's principle that 

H £ (t,0 = e^ (0 - i fe^ [ m ~ ^ dfrUft - v,t)dv 

Jo JR d ?1 

:= e^ m uo(0 + A E Lle(t, £) + E £ U £ (t, £), 
where the operator A £ is defined in (4.1) and may be recast as 

t V 

A £ ii £ = - [ e^ v [ e ls ^- em ^d Sl R(Ci - eZ)d£iUe(t - v, $dv 

Jo Jm. d Jo 

Jo i(£?-£ m z m y 

and the remainder E £ is then given by 



(i^+e- P e(o)ut,o = o (42) 
ii e (o,e) = u (e), 



(4.4) 



(4.5) 



1 



E ^e =i / ^ e ^^i " ^)^iH £ (i - v, £)dv 

Jo JR d ?1 ~~ e ? 

+ * 1 1 rf^y (e "" " e ^ )i?(6 " emMt ~ v > (4 ' 6) 
■=h+h. 

16 



(4. 



For the calculation of I\, note that equation (4.2) has the explicit solution: 

U(t,0 = j^-KWuoit). (4.7) 

On one hand, we may obtain the following expression of integral in v using the method of 
separation of variables: 

JO si ^ s JO 

Jo 

On the other hand, we have the following simple estimate of this integral: 

[* e ^(tT-c m Z m )iX £ (t - v,£)dv < f le^f-^iX^t - v,0\do < C\u (O\- (4. 
Jo Jo 

This gives that 
Jo 

We therefore find the estimate of using Lemma 5.9: 

\h\<Ce m ft 1 A 1) maxdfloffll, If flo(fl|}£(ft ~ 

jR d Si _ e ? e (4.11) 

< Cmax{l,e m }max{e m ,^- m | loge|(f + l)i~ 2 , e d ~ m (f + 1) ^~ 2 } |^ (0 1 - 



< ^( m \ m ,m A ^)max{|uo(C)|, If «o(0l}- (4-10) 
si £ s 6 



For the calculation of I2, we first estimate the integral in v. 
e m / -7s -— H £ (t-v,0^ 



< ^( m Lfm A ™«{l«0(0|, If *>(0I}- ( 4 - 12 ) 
Si 6 s e 



Hence 



1*1 < C £ff(fr^ (413) 
<Cmax{l,nmax{ £ m ,^- m |log £ |(r + l)-^,e d - m (r + l)-" 2 }irno(6l- 

Note that f 

\A £ U(t,0\<C [ \U(s,0\ds, (4.14) 
J 

over a bounded interval in time. The equation 

(I-A £ )U(t,0 = S(t,0 (4.15) 
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therefore admits a unique solution by Gronwall's Lemma, which is bounded by 

\U(t,0\ < ||S , || 0O e c *. (4.16) 

We verify that the solution to 

(I-B £ )H £ = e u ^u (O, (4-17) 

is given by 

H e (t,O = e^ m -^u (£). (4.18) 
The error V e (t,£) = (f7 £)S (t, £) — il £ (i,£)) is a solution to 

(I-A e )V e = E e ii e (t,t), (4.19) 

so that over bounded intervals in time, we find that 

\U EtS (t,0-He(t,0\ = \V £ (t,0\ <max{ £ m ,^- m |log £ |}. (4.20) 

□ 

From our assumption that R(£) 6 C 2 (M. d ), we find that 

| e «(e»-M0) _ e it(r- P )\ < t \ PE ^ _ p) \ < ^2 (421) 

The reason for the second-order accuracy is that Jf2(— £) = and Vi?(0) = so that 

first-order terms in the Taylor expansion vanish. 

In terms of the solutions of PDE we defined in (1.3), we may recast the above result as 

\(U £yS -U)(t)\ < max{eW rf - m | loge|} W(t,0 = e^"^^), (4.22) 
We now prove Theorem 1.1. By the triangle inequality, we have the estimate 
E\\u £ (t,0-U(t,0h 

no(e)— 1 +oo +oo 

<E|| J2 («.,„ - Ul s ){t)\\ 2 + E|| ^ n £ , n (t)|| 2 + \\{U £tS -U){t)\\ 2 + || ^ U- S (t)h- 

n=0 no(e) n (e) 

(4.23) 

The vanishing of the first three terms on right hand side of this inequality when e goes 
to zero follows from (2.53), (3.35) and (4.22) respectively. The fourth term also vanishes 
because of the L 2 convergence of U s a . 



5 Inequalities and Proofs 

In this section, we present and prove several inequalities used in earlier sections. There are 
similar versions of Lemma 5.2 and 5.7 in [6]. The proofs are given below for the convenience 
of the reader. The proofs of similar versions of Lemma 5.1 and 5.8 can be found in [6]. 
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Lemma 5.1. We have the following identity for r/ > 0: 

We also claim the following estimate with n := \I\ — 1: 

t" 



+ ir) 



(5.1) 



\K(t,£,I)\< 



in 



Lemma 5.2. Assume rj > 0. We /iaue i/ie following inequality: 



da 



< 



C 



\a + A + ir]\\a + B + in\ ~ \A - B + in\ 
where log + x := max{0, log x} for x > and log + := 0. 



1 + log. 



A-B 



(5.2) 



(5.3) 



Proof. Without loss of generality, we assume B > A. We split the integration over 
(— oo, — A \ B ) as follows: 



A+B 
2 



A+B 
2 



da 



\a + A + in\\a + B + vq\ 
da 



+ 



da 



+ 



B \a + A + ii]\\a + B + ir/\ J - B -(B-A) \ a + A + irj\\a + B + ir]\ 

~B-{B-A) da 



\a + A + iry||a + 5 + ^| 
The first term is estimated as 

da 



A+B 
2 



< 



B 



a + A + irj\\a + B + irj\ \ - + A + irj\ J |a + -B + ^ 



A+B 
2 



da 



< 



B-A 

2 / 2,, rfa 



|A - B + ir/| 7o Va^TT 



< 



C 



\A - B + zr/| 



1 + log H 



A-B 



Likewise, the second term is estimated as 



da 



< 



B~{B-A) \a + A + in\\a + B + in\ ~ \A - B + 



-B 



da 



< 



1 



< 



\A- B + in| 
C 

| A - 5 + Z77| 



ja + S + ir/l 

a-A 
i da 



1 + log_ 



Va^Tl 
A-B 



V 



(5.4) 



(5.5) 



(5.6) 
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We obtain the bound for the third term by using the inequality \a + A + vq\ < \a + B + irj\ 
on (-00, -B - (B-A)). 



-B-(B-A) 



d,a 



la + A + iri\\a. + B + 



< 



-B-(B-A) 



da 



00 da 



r] J b-a a 2 + 1 



If I? — A > r], we have 



+00 



da 



b-a a 



2 + 1 



< 



x da r] 

3_A ^2 - 5 _ A' 

>7 



in which case 



-B-(B-A) 



If B — A < r], we have 



da 



< 



1 



< 



V2 



\a + A + irj\\a + B + in\ ~ B - A ~ \A - B + irj\ 



b-a a 2 + 1 



< 



da 



+ 



G?a 7T 



a 2 + 1 I b-a a 2 + 1 



<T + (1 



5 - A, 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



in which case 



-B-(B-A) 



da 



I, it . \/2(7r/4 + l) 

< -( T + 1) < 



|a + A + ir/||a + 23 + irj\ '' 77M ' ~' ' \A-B + ir]\' 
By symmetry, the integration over (—00, — A \ B ) admits an identical bound. 



(5.11) 



□ 



Proposition 5.3. Assume r],fj > 0, fj > 2rj, and n 1 bounded. We have the following 
inequality: 



1 



1 



C 



\a — a + i(n — fj)\ |d + A + in\ \a + B + zf/| 



dada < 



1 + log. 



A—B 



\A - B + Mjl 



(5.12) 



Proof. Without loss of generality, we assume 73 > A. 
We first integrate in a by using Lemma 5.2: 



1 



da < 



1 



d — a + z(r/ — f/)| |d + ^4 + irj\ ~ \a + A + ifj\ 



1 + log+ 



a + A 



'1 



(5.13) 



The integration over (—00, 



A+B^ 



is split as in (5.4) into three pieces: (— B, 



(B - A), -B) and (-00, -B - (B - A)). 



A+B \ 
2 h 



-B 



The term 1 + log_ 



is bounded by 1 + log. 



integral over (— B, — A ~t B ) is then bounded by 



A—B 



n 



on the interval (-B, - A ^ B -). The 



A+B 
2 



1 +10S 



a+A 



da < 



C 



\a + A + ifj\\a + B + ifj\ \A — B + ifj 



1 + log. 



A- B 



(5.14) 
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The same bound holds for the integral over (-B — (B — A), —B). 

It remains to estimate the integral over (—oo,—B — {B — A)). The domain of integration 



a+A 



can be divided into two parts. On the set \a + A\ < r), we have log_| 

applying Lemma 5.2 immediately gives the desired result. On the set \a + A\ > fj, we use 
the simple bound on this interval 



0, so that 



1 + log 

It then follows that 

-B-(B-A) 1 + log 



a + A 



< 



a + A 



V 



+ i 



1/2 



(5.15) 



a+A 
V 



-^da < C 
= C 



-B-{B-A) 



oo 
+oo 



\a + B + ifj\ 3 / 2 



da 



da 



\a + A + ifj\ \a + B + ifj 

-r r 

-a (a 2 + 1)3/2' 

After performing an analysis similar to that from (5.8) to (5.11) we find the estimate: 



(5.16) 



-B-(B-A) 1 + log. 



a+A 



7 da < 



C 



1 + log , 



A-B 



V 



\a + A + ifj\\a + B + ifj\ " \A - B + ifj\ 
which concludes the proof. 

Lemma 5.4. Assume rj > and r/ _1 is bounded. We have the following inequality: 

d£ < Ce x (a^ V 1). 



(£) 2d \e m a + i m + ie m n\ 2 



(5.17) 
□ 

(5.18) 



Proof. The domain of integration may be split into two parts. On the set \e m a + £ m | > 1, 
the integral is bounded by e m ||l/(£) 2rf ||i. On the set \e m a + £ m | < 1, we change to polar 
coordinates 



-di<C 



£ 



d-l 



| e m a+e »|<r (0 M |e m a + e , + te m »/ 12 



d\t\. 



l\e*«+e*\<\ ~ (0 2d k m a + C m + ie™ V \ 2 
When m < d < 2m, we use the inequality \l/(^) 2d \ < 1 and change variables to Q = e m a+£ n 

'-d-l 
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L 



hgm|<r \e m a + £ m + ie m r]\ 2 



\Q + ie m r)\ 2 



(5.19) 



(5.20) 



On the set [0,£ m ?7], we bound \Q + ie m r]\ from below by its imaginary part and estimate the 
integral 



" m (Q-e m a 



-dQ < 



'"VP*- 1 + ^-\ Q < ce^ja V (5.21) 



\e m r]\ 
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On the set [e m ry, 1], we bound |(5+ie m f?| from below by its real part and estimate the integral 







Q 2 







Q- 



-dQ <Ce d - m {aVl)™-\ (5.22) 



This concludes the proof for m < d < 2m. When d > 2tn, we use the inequality \£ d 2m / < 
1 instead of |l/(0 2rf | < 1 and the rest of the proof is similar. □ 

Proposition 5.5. Assume i] > and ri~ l is bounded. We have the following inequality: 

„ 1 1 



(0 2d \e m a + £ m + ie m r]\ 2 



d^Ce^^a^Vl). 



(5.23) 



Proof. To prove (5.23), we just need to use l/(£) 2d to bound an additional ^( d ~ m ) s and the 
rest of the proof is the same as in Lemma 5.4. As a result, we lose a contribution e( rf ~ m ) 5 
for lowering the order of a by □ 



Proposition 5.6. Assume n > 0, ?] 1 bounded, and k £ N. VFe Ziaue £/ie following inequal- 
ity: 



1 



1 + log H 



-d£ < Ce x {a V 1) 



(5.24) 



Proof. The integral on the set 



a+g m /e" 



Lemma 5.4. The integral on the set 
Lemma 5.7. Assume ?/ > 0. We have the following inequalities: 

1 1 



< 1 is immediately bounded by Ce am using 



> 1 is bounded by Ce m ||l/(£) 2d lli • 



sup 



\a + C m + iri\(Z - uj) 2d {C) 2d ~» ^ ~ {a} ' 
where the (a) in the denominator again has to be dropped if we remove the factor (£,} 2d 

Proof. We will first show that 

1 1 



(5.25) 



|a + £ m + n?l {i-uj) 2d 
If rj > 1, the integral is bounded by 



d£ < C\\ogr]\. 



1 



\a + £ m + irj\ (£ - uj) 



2d 



d£ < 



1 



l —dt<-<C\\ogr 1 \. 



V 



(5.26) 



(5.27) 



Therefore, we only need to look at the case for rj < 1. The integral over \a + £ m | > 1 is 
bounded by 



i 



i 



a+ £m|>l \a + i m + irj\ 



2(/ 



d£ < C 



1 



2,1 



d\i\ < C < C\\ogr]\. (5.28) 
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The integral over |a+£ m | < 1 can be estimated by splitting the integration domain according 
to the size k<\^ — uj\<k + l for k = 0, 1, • • • 



1 ^ 



K m|<i \a + i m + ir]\ 

: ^ /|a+e»|<i}rHfc<|€-o;|<fc+i} + e + {t-u)™^ (5-29) 



=E / yi / J(ia ; ! 1 --; 2 f- i) ^-^-, 

J\a+?«\<1 | a + C m + «?| ./fc<|£-a,|<fc+l (£ ~ w ) 

where J(|£|, #i, • • • j^d-i) denotes the Jacobian to change to polar coordinates. For fixed 
|£|, we have the inequality 

/ 7t \2d ' ' ' " y d-l - / /-, , j 2 \d dUl ' ' ' d& d-l 



(i+/c 2 ) d 



c^ier 1 , lei < i. 



Summing up the terms over k = 0, 1, • • • gives 



(5.30) 

<(7|log77|. (5.31) 



We now prove (5.25). Without loss of generality, we may assume \a\ > 2. 
The integral over the domain \a + £ m | > a/2 is easily bounded by 



+^I<M |q + e m + ir/| <£ - (0 M " l«l J\c+?*\<&. (e - (0 

< g < C| log 7/1 
~~ |a| ~~ (a) 



(5.32) 



On the domain \a + £ m | < \a\/2, we have £ m > |a|/2. Note also that (£) 2d > (£ m )- The 
integral over this domain is therefore bounded by 



L+^\<M la + e m + «7| <£ - w) M (0 M - (a) L 



1 « 



C|logr/| 



< 



(a) 



(5.33) 
□ 



Lemma 5.8. Let 

(Y-,m(n ■.= f - 



■■= I ^"ST^ (5-34) 
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be a family of linear operators parametrized by a complex parameter z = a + in with n > 0. 

(i) For al > 3 we have 

\Y Z R\ < C\\R\\ 2d>2d . (5.35) 

(ii) If z' = a' + irj and i] > rf , then for d > 3 

\Y Z R - Y Z ,R\ < C\z - z\r,-^-^\\R\\ 2dM . (5.36) 

Lemma 5.9. Let R satisfy the smoothness condition defined in Section 1. We have the 
following inequality: 

e ' n / ( ^^y A ^^r^ i < Cmax^.e^llogeKr+lJ^^Cr+l)-- 1 }. 

(5.37) 

Proof. We decompose the integral into three integrals on > e m ^ m + 1}, {e m (^ m + 1) < 
< £ m^m + jj and {o < £™ < e m (£ m + 1)}. Clearly, the first integral in bounded by 



to 



£ m j _L«_£i2^ 1 _ p or seC ond integral, note that R is bounded, and use ^m_*mgn 
bound the term in the bracket on the left hand side of (5.37). We then change variable to 
polar coordinates and let Q = £™ to obtain 



e- f ( 1 A V^-^V < Ce- ! Qm " dC 

<Ce d ~ m \loge\(C + l)i~ 2 . 

(5.38) 

For the third part, we use l/e m to bound the term in the bracket on the left hand side of 
(5.37) and apply the same type of change of variable to obtain 



e m / (- — A — ) — — —d£i <t Q™~ z dQ 

'{o<?r<£ m (C m +i)} ?i e ? e si Jo (5.39) 

< C£ d - W (e m + 1)™ -1 . 

□ 
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